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ABSTRACT 

Context. The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster 

and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. 

Aims. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the difl"erent heating mechanisms. 

Methods. HIFI observations of mid-/ transitions of CO and HCO^ isotopes allow us to bridge the gap in excitation energies between observations 

from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed 

line profiles allows to identify the physical structure of the different components. 

Results. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to ^ 2 THz can 

be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. 

We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of 

cooler, but nevertheless dense, molecular clumps. 

Key words. ISM: structure - ISM: kinematics and dynamics - ISM: molecules - HII regions - Submillimeter 



1. Introduction 

DR2 1 is a deeply embed ded H n region created b y the radiation 
from at least six OB stars (iRoelfsema et aliri989h . It sits within a 
ridge of dense molecular material that obscures the H ii-region at 
optical wavelengths. The embedded cluster drives a violent bipo- 
lar outflow in north-east to south-west direction. It is prominent 
in the 2 fim emission of vibrationally excited H 2, tracing hot , 
shocked gas and i n the Spitzer 4.5 fim channel JGarden et al.L 
Il986t iDavis et all 12007 !). As the cluster is located close to the 
eastern edge of the molecular ridge, the eastern, blue-shifted out- 
flow expands in a blister-like fount ain, while the we stern, red- 
shifted outflow is highly coUimated (jLane etal.Lll990 ). 

Spitzer 8 jum images r eveal spots of brigh t PAH emis- 
sion with a size below 10" ('M arston et all 12004) . They repre- 
sent the surfaces of high-density, UV irradiated clumps, form- 
ing photon-dominated (or photo-dissociation) regions (PDRs), 
transition zones from ionized and atomic gas to dense molec- 
ular gas where physics and c hemistry are domina t ed by 
UV radiation from young stars jHoUenbach & Tielensl 119991: 
lOsse nkopf et al., 2007). These hot and dense regions give rise to 
the emission of PDR tracers, such as HCO^, high-7 CO, atomic 
and ionized carbon, and atomic oxygen. 

An overvi ew on the numerous e xistin g observatio n s in D R21 
is given by Schneider etaP (120061) and IJakob et al.l (l2007h . In 
spite of the wealth of data, the heating of the molecular gas is still 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



debated. The bright H2 emission indicates shock heating in the 
outflow and the wings of the molecular lines prove its dynam- 
ical impact, but excitation m odels of the observed emission of 
CO, [Ci], [Cii], and [Oi] bv lLaneeta l. (1990) and LTakob et all 
(l2007h have shown that the emission of those tracers cannot be 
explained by shocks, but is consistent with a pure UV heating, 
i.e., PDR physics. 

To quantify the heating of the gas, we analyze spectra taken 
with the HIFI instrument (de Gr aauw et al., 2010) o n board 
the Herschel Space Observatory dPilbratt et al.L 1201 Ol) during 
the performance verification campaign. Observations of high-/ 
HCO^ t ransitions trace hot material, i onized by UV radiation or 
X-rays (ISternberg & Dalgarno'. 1995). Hot water lines are pro- 
duced in shocked gas (Snelletal., 2005), and observations of 
CO isotopes around 7=10 close the gap in the excitation lad- 
der between ground-based observations and existing ISO data, 
allowing to obtain a full picture of the temperature distribution. 

In Sect. 2 we present the observational data. Sect. 3 com- 
pares the measured line profiles to distinguish different compo- 
nents based on their velocity distribution. In Sect. 4 we provide 
a model for the emission, supporting the PDR character of the 
source, and discuss the results in Sect. 5. 



2. Observations 

2.1. HIFI observations 

All spectra presented here were obtained in performance verifi- 
cation observations for the HIFI instrument. As their main goal 
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Table 1. Summary of the used HIFI observational data 



transition 
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HCO+ 6-5 


535.062 
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OTF map 
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0.04 


HCO+ 12-11 


1069.694 
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FSW spectral scan 
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'^CO 10-9 


1101.350 


21 


LC spectral scan 
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0.08 


C"*0 10-9 


1097.163 


21 


LC spectral scan 


270 


0.08 


[Cii] 


1900.537 
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DBS raster map 


14 


1 






20 


convolved-l-binned 


112 


0.2 






40 


convolved-l-binned 


420 


0.1 



' DBS = dual-beam-switch, OTF = On-The-Fly, FSW = frequency- 
switch, LC = load-chop, OFF position = 20h37ml0s, 42°37'00" 



was to demonstrate the functionality and performance of the dif- 
ferent observing modes, the spectra were taken with a large vari- 
ety of observing modes and strategies. Consequently, every spec- 
trum was taken in a slightly different manner. All observational 
parameters are summarized in Table [T] 

Most observations where single-point observations to- 
wards the central position of the DR21 Hii region at 
RA=20h39m01.1s, DEC=42°19'43.0" (J2000). Fully sampled 
maps were only obtained in the [CII] line. Data were taken with 
the wideband spectrometer (WBS) at a resolution of 1 . 1 MHz, 
corresponding to 0.2 km/s (at 1900 GHz) - 0.7 km/s (at 
500 GHz). The [CII] data were rebinned to a velocity resolution 
of 0.45 km/s to improve the signal to noise. 

2.2. Complementary data 

ISO Long Wavelength Spectrometer 43 - 197 yum grating scans 
were obtained for the DR21 central position from the ISO Data 
Archive (TDT 15200786). Integrated line intensities were ex- 
tracted for [Oi] at 63 and 145 yum and the CO 14-13 to 17-16 
transitions. Mid-7 CO lines of the DR2 1 region were ma pped 
with the KOSMA 3 m submm telescope dJakobet alil2007n . We 
use the lines of CO and '^CO from J = 3-2 to 7-6, which have 
been observed at native angular resolutions from 80" to 40". The 
HCO+ 1-0, Hi3CO+ 1-0, and HC O+ 3-2 observations we re taken 
with the IRAM 30 m telescope ( (Schneider et al.L l2010h . Native 
angular resolutions at the 1-0 and 3-2 transition frequencies are 
28" and 9", respectively. 

2.3. Beam size effects 

For a direct comparison of the different data sets we smoothed 
the available data to the coarsest common angular resolution of 
80", matching that of the CO 3-2 KOSMA beam and the lowest 
frequency ISO observations. This is impossible for the single- 
point HIFI observations. Moreover, the Herschel beam varies 
between 40" HPBW for the HCO+ line at 535 GHz and 20" for 
'^CO at 1100 GHz. Those data were corrected for the different 
beam filling by estim ating the sou rce si ze from th e HCO^ 3-2 
line as a PDR tracer dSternberg & Dalgarnoi Il995i) with quite 
compact emission (20" to 30" in diameter). Successively con- 
volving from the native angular resolution of 9" to the angular 
resolution of the Herschel data, and further to the final, smoothed 
spectra at 80", we obtained scaling factors for the HIFI spectra, 
being 0.5 when going from 40" to 80" and 0.33 when going from 
20" to 80". The maps obtained in [CII] allow a direct smoothing 
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Fig. 1. HIFI spectra of the lines used in the fit of the PDR prop- 
erties. The [Cii] data are convolved either to the resolution of 
the HCO+ 6-5 line or to that of the CO 10-9 lines. The negative 
features in the HCO^ 12-11 line are artifacts from the frequency- 
switch observing mode. 



to a resolution of 40" (see Table[TJ. Beyond that size, the same 
scaling factor as above was applied. 



3. Line profiles 

The molecular ridge including the DR21 Hii region has an in- 
trinsic LSR velocity of -3 km s"' . The qu iescent material is vis- 
ible in narr ow absorption lines of NH3 (iMatsakis et al.L Il977h 
and H2CO ("Biegin g et al.L fl982h . A second velocity component 
at 8-10 km s"' is known to be associated with the W75N com- 
plex. It appears in emission in CO and [C i] (Jakob e t an,l20()7h . 
as a narrow absorption feature in HCO^ 1-0 (.Nvmanill983 l). and 
has a very broad velo c ity distribution in the Hi 2 1 cm absorp- 
tion dThompson et al.L 119691: iRoberts et all Il997h . The wings 
of the low-/ CO lines trace outflow velocities down to about 
-20 km s"' for the eastern, blister outflow and up to w 10 km s"' 
for the western outflow. 

Figure [T] shows the profiles of the HIFI spectra of CO iso- 
topes, HCO^ and [Cii]. All lines peak at about -4 km s '. The 
CO and HCO^ lines have a similar shape, but the [Cii] line 
shows an additional broad blue wing extending down to -30 km 
s"'. This indicates that the warm molecular material is slightly 
blue-shifted relative to the cold gas and that the [C 11] emission 
is not only originating from that warm gas, but also from the 
ionized wind in the blister outflow. 

Figure |2]compares the shapes of the CO and HCO^ isotope 
lines with complementary ground-based measurements towards 
the same positions. We show only a few selected transitions as, 
e.g., the data for the CO 7-6 or 4-3 lines provide no additional 
information. All CO isotopic lines up to 7-6 are roughly symmet- 
ric, centered at the ridge velocity of -3 km s '. The lines of the 
main isotope are heavily self-absorbed with the absorption dip 
marking line center. The 10-9 lines, tracing hotter material, are 
slightly asymmetric and shifted to -4 km s '. All HCO^ lines, 
except H'^CO^ 1-0, have a profile very similar to the 10-9 lines 
of the CO isotopes indicating that HCO^ is mainly abundant in 
the heated layer and less in the overall molecular material. The 
H'^CO^ 1-0 line has a deviating profile centered at -2 km s ' 
tracing global infaU (Kirby, 2009). CO 3-2 and HCO+ 1-0 show 
an additional feature at the W75N complex velocity. 
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Fig. 2. Selected profiles of CO isotope lines (upper plot) and 
HCO"^ isotope lines towards the DR21 central position. All lines 
are convolved or scaled to the coarsest common resolution (see 
Sect|231l. 



To better understand the exact velocity distribution in the 
source, we pl ot some additional HIF I hnes towards th e same 
position (see Falparorie et aU 120101 : lOssenkopf et all l2010t 
Ivan der Tak et all |2010 S) in Fig. [3] The 1097 GHz line of hot 
water matches the rest of the PDR tracers. The foreground at - 
3 km s"' is visible in absorption in the ground state transitions 
of water, but in emission in a nearby OFF position for [Cii]. 
At velocities of about +2 km s"', we find the warm gas facing 
the redshifted western outflow, apparent as a secondary peak in 
the ground state water and CH^ lines and as a shoulder in '''CO 
and the hot water line. Finally, we can clearly identify the cold 
W75N component at ^ -9 km s"' . 

For the warm gas we can distinguish two velocity com- 
ponents - a blue shifted one related to the blister outflow at 
-4 km s"' and a second one, at 2 km s"', related to the west- 
ern outflow. This is consistent with the clumpy PDR geometry 
proposed bv'Lane et al.l (1199 0'). In terms of line modelling, it is, 
however, impossible to separate the two outflow directions as 
their emission is overlapping to a large degree. We simply add 
their intensities. The [C ii] emission requires a special treatment 
because of the additional emission from the ionized outflow. To 
take this into account, we have obtained an integrated [C ii] in- 
tensity by scaling the '^CO 10-9 line profile to match the peak 
and the red wing of the [C ii] profile and ignoring the remaining 
blue-wing emission. 
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Fig. 3. Auxiliary lines measured by HIFI towards the same posi- 
tion. They can be used to estimate the contribution of the fore- 
ground material. The '^CO 10-9 profile is also displayed for 
comparison. 



4. Modelling 



We use the KOSMA-t PDR code dRolhg et all l2006h to model 
the emission of PDR ensembles, representing a distribution of 
spherical clumps with dN/dM oc M"'** (Cubicket al., 2008|). 
For DR2 1 two ensembles with different properties had to be su- 
perimposed, a hot component, close to the inner H ii region with 
strong FUV illumination, but only a small fraction of the to- 
tal mass, and a cooler component that fills a larger solid angle 
and provides the bulk of the material. Each clumpy PDR en- 
semble has five free parameters: the average ensemble density, 
Hens, the ensemble mass. Mens, the UV field strength, x given in 
units of the Draine field, and the minimum and maximum mass 
of the clump ensemble, [M^in, M^ax]- In contrast to most other 
PDR models, we fit absolute line intensities, using the available 
ground-based observations, complementary ISO data, and the 
HIFI lines of the CO isotopes, HCO^, atomic and ionized car- 
bon, and atomic oxygen. The chemical network that has been 
applied in these calculations includes '^C but not ''^O. The C'^O 
lines were scaled from the '^CO intensities with a conversion of 
1:8. Simulated annealing was used to find the optimum parame- 
ter combination. 

The significance of the model is limited by the fact that the 
clump superposition ignores mutual line shading between differ- 
ent clumps, i.e., optical depth effects are only considered within 
individual clumps. This is usually justified by the virialised ve- 
locity dispersion between different clumps, but for optically very 
thick and broad lines some correction is needed. To estimate the 
effect we have computed the optical depth for the bulk of the 
individual clumps. This is of the order of unity or below for the 
majority of the observed transitions, reaches values up to ten for 
'^CO and HCO^ transitions up to 7 = 4, but exceeds ten for the 
CO main isotope lines up to 7 = 6 and the [O i] 63 fj.m line. For 
the CO lines showing clear self-absorption dips, we performed 
a Gaussian fit to the line wings and used the integrated inten- 
sity of that Gaussian to compute the total emission including the 
blocked radiation from the inner clumps close to the H ii-region. 
As we have no spectral information for the [O i] line, we have 
no estimate for the blocked radiation in this case, so that we ex- 
cluded that data point from the fit. 

The best fit result is shown in Figure H) The corresponding 
model parameters are: 
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Fig. 4. Two-ensemble PDR model fit to the observed CO, HCO+, 
and fine structure line intensities, shown as function of the upper 
level energy. HIFI measurements are depicted as open triangles, 
complementary data points as open circles. The dotted line and 
the shaded areas indicate the contributions from the two ensem- 
bles to the CO and the fine structure lines, respectively. 

Ensemble 1: ;t' = 1.0 10^ «ens = 1-3 10''cm-\ Mens = 150 Mq, 

[M„i„,M„,ax] = [10-2,8 10']Mo 
Ensemble 1: x ^ 3.2 10^ Wens = 1-1 10^cm-\ Mens = 830 Mq, 

[M„i„,M„ax] = [10-3,10i]Mo 
The two-ensemble model provides a reasonable fit to all 22 data 
points. The model slightly overestimates the emission for the 
optically thick rotational lines of '-^CO and HCO^ which showed 
no clear self-absorption dip, so that we have insufficient optical 
depth corrections. The reduced x^ value of the fit amounts to 
11.1, where a contribution of 9.5 results from the three optically 
thick low-y data points of '-^CO and HCO^. For the remaining 19 
data points we would obtain an excellent reduced;^'^ value of 1 .6. 
The topology of the x^ function shows several side minima, but 
they are worse by at least a factor two. As expected, the model 
predicts a too high [Oi] intensity, as it ignores that the outer 
clumps of the cooler ensemble block the contribution from the 
hot inner component. 

5. Discussion 

The double-peak structure of the predicted CO intensities across 
the excitation ladder reflects the two different UV fields leading 
to different excitation conditions. No single-parameter ensemble 
can fit all data, but distinct low temperature and high tempera- 
ture components are needed. The parameters of the model are in 
agreement with independent estimates. A UV flux of 10^ Draine 



fields corresponds to a geometrical distance of 0.06 pc, i.e. 7", 
from the central cluster, matching the size of the PAH emis- 
sion (Marston et al. , 2004). Our clump densities match those de- 
termined by Jako b et al.l (l2007h for the extended cool gas, but 
are sligh tly higher than the i r hot-gas density ( 4 10^ cm'^). In 
contrast, iJones et alj (11994 : Ivan der Tak et alj (1201 Ol) find still 
somewhat higher densities for the hot gas, up to 10^ cm"^. The 
hot ensemble mass is cl ose to the 170 Mp derived from early 
CO 7-6 observations by iJaffe et all (Il989l) . The total mass of 
the PDR ensembles falls between the mass limits derived by 
I Jakob et a l. (2007) from dust observations and from line radia- 
tive transfer fits. 

While the existing ground-based observations provide a very 
good constraint on the properties of the extended cool gas, and 
the ISO lines show the total amount of hot gas, it is only the 
set of new HIFI data that puts the hot and cold distributions 
well a part from each othe r in terms of the temperature structure. 
While lJakob et alj (l2007h obtained cooling curves with single 
peaks, the new data for the 10-9 lines of the CO isotopes and the 
HCO^ transitions force the fit to a bimodal distribution of exci- 
tation conditions. When we exclude the Herschel data from the 
model fit, we obtain a parameter set that shows a UV field that 
is lower by a factor ten for the hot ensemble, i.e., that would im- 
ply molecular clumps farther away from the central cluster. For 
the cold ensemble, the fitted UV field is also somewhat lower, 
while all other parameters remain similar to those from the full 
fit. Only with the Herschel data, we therefore obtain a parameter 
set that is consistent with the source geometry. 

As the two-ensemble PDR model is able to fit all of the ob- 
served lines, we find no evidence for a shock heating of the dens e 
gas. This is in agreement with the analysis of lLane et alJ(ll990l) . 
explicitely excluding a shock origin of the fine-structure Unes, 
but seems to be in contradiction with the analysis of the line pro- 
files in Sect.[3]that shows excited outflow material. We conclude 
that the material visible in the blue line wing, characterizing the 
blister outflow, is contained in dense clumps that are accelerated 
by the outflow, but that are chemically and energetically fully 
dominated by the UV field and not by the associated shock. 
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